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We derive a simple empirical photometric redshift estimator using 
a training set of galaxies with multiband photometry and measured red- 
shifts in the Hubble Deep Field (HDF). This estimator is model- independent; 
it does not use spectral templates. The dispersion between the estimated 
redshifts and the spectroscopically measured ones is small; the disper- 
sions range from cr^ ~ 0.03 to 0.1 for z ^ 2 galaxies, and from o"^ ~ 0.14 
to 0.25 for z ^ 2 galaxies. The predictions provided by our empirical red- 
shift estimator agree well with recently measured galaxy redshifts. We 
illustrate how our empirical redshift estimator can be modified to include 

\J^ ' flat spectrum galaxies with 1.4 ^ z ^ 2. 
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1. Introduction 



Our aim is to derive a simple but accurate empirical photometric redshift esti- 
mator which is mo del- independent. We do not use spectral templates; we only 

O I use a training set of galaxies with multiband photometry and measured red- 

shifts. We have applied our method to the Hubble Deep Field (HDF) (Williams 

jrt ! et al. 1996) and derived a catalog of estimated redshifts for 848 HDF galaxies 

with I < 27 and with measured fluxes in U, B, V, and / (Wang, Bahcall, & 
Turner 1998). 
^ ' To derive formulae for the estimated redshifts, we first divide the galaxy 

redshift sample (the training set) into regions of high and low redshifts (z '^ 2 
and z ^2) based on empirical color cuts. We then divide both regions into color 
ranges to minimize dispersion. The physical motivation for this latter step is 
to reflect color shifts with z for different type galaxies. We then find the best 
linear fit between redshift and colors for each color range, 

Za = Cl+ C2{U -B)+ C3{B -V)+ C4(y - /), (1) 

where Cj {i = 1,2,3,4) are constants. Here we have assumed that photometry 
has been done in four bands, U, B, V, and /. The dispersion (a^) between our 
estimated redshift, Za, and the measured spectroscopic redshift, z, is calculated 
using the jack knife method. For a sample with A^ data points, we make A^ 
subsamples, each omitting one data point. We then carry out linear fits A^ 
times on the A^ — 1 data point subsamples. We use each fit to make a prediction 
for the omitted data point; the rms of the true values of the omitted points 
about the fit predictions, scaled by [(A^ — 1)/A^]^'^, is the rms dispersion of our 
fit. The estimated dispersions of our formulae are therefore quite robust. 
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The error in Za due to photometric errors is 

AZa 



ic2AUf + 1(C3 - C2)A5[2 + |(C4 - C3)AV\'^ + (cAAlfV^^ . (2) 



2. Application to the HDF 

We used 82 galaxies with measured spectroscopic redshifts z ~ 0.1 — 3.5 in the 
HDF as our training set (Cohen et al. 1996; Hogg et al. 1998; Lowenthal et al. 
1997; Phillips et al 1997; Steidel et al. 1996). We did not use all 90 (excluding 
three z > 2 galaxies with erroneous/uncertain z's) HDF galaxies with measured 
z's. The eight outlying galaxies have az = 0.45, mostly near the boundaries of 
the color ranges; these galaxies were not used in deriving our redshift estimator. 
We found that for z '^ 2, the galaxies satisfy at least one of the following 
three color selection criteria: 

U > 25.66, U -B> 0.91, B-V < 1.37, V -I < 0.5; (3) 

I > 23.5, U - B > 2.2; (4) 

I >23.5, B-V >2.2,U -B > -0.5. (5) 

The galaxies with z < 2 in the training set do not satisfy any of the above rela- 
tions, suggesting that z < 2 galaxies generally fall outside these color- magnitude 
regions. 

We divide z < 2 galaxies into three color ranges (cr), and determine the 
best-fit redshift estimate formula and dispersion for each color range: 

1. cr=l: {U -B) <{B-V)- 0.1 (28 galaxies) 

Za = 0.4111-0.1852(C/-S)-0.3062(S-y)+0.7301(y-I), a^ = 0.034; (6) 

2. cr=2: {U - B) > {B -V) - 0.1 > {V - I) (21 galaxies) 

Za = 0.163 -0.171{U-B) + 0.3iO{B-V) + 0.194:{V- 1), a^ = 0.095; (7) 

3. cr=3: {U -B)>{B-V)- 0.1 < (V - I) (19 galaxies) 

Za = 1. 126 + 0.480{U-B)-0.513{B-V)-0.250{V- 1), cr^ = 0.097. (8) 

We divide z >2 galaxies into two color ranges, and determine their redshift 
estimate formulas and dispersions: 

1. cr=4: (B-V)- 0.5 > {V - I) (8 galaxies; z^3) 

Za = 2.37 + 0.02{U -B) + 1.61{B -V)-2.47{V -I), a^ = 0.14; (9) 

2. cr=5: (B-V)- 0.5 < (V - I) (6 galaxies; z ^ 3) 

Za = 2.18 + 0.10{U -B) + 0.20{B -V) + 0.75{V -I), a^ = 0.25 (10) 

Fig.l shows our estimated redshift Za (given by Eqs.(|6|)-(p^) versus the 
spectroscopic redshift z for 90 HDF galaxies. 
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Figure 1. Our empirically estimated redshift Za (given by Eqs.@- 
(ffH)) versus the spectroscopic redshift z for 90 HDF galaxies. The 
galaxies with known measurement errors in UBVI are plotted with 
error bars in Za- Error bars without points represent the eight galaxies 
not used in fitting. 





Figure 2. (a) Our empirically estimated redshift Za (given by Eqs.@- 
(0)) versus the Sawicki et. al. (1997) template-fitting photometric 
redshift ztemp, for 848 galaxies in the HDF with I < 27 and measured 
UBVI. The symbols are the same as in Fig.l. The solid diagonal line 
indicates Za = ztemp] the dotted lines mark the region \za — ztemp\ < 0.5. 
(b) Same as Fig.2(a), but with the addition of Eq.([TT|) to select flat 
spectrum galaxies with estimated redshifts given by Eq. (|l2|) . 



Our empirical photometric redshift estimator is consistent with template- 
based photometric redshift estimators where calibrating redshift samples are 
available. Fig. 2(a) compares our empirically estimated redshift Za (given by 
Eqs.(P)-(|lO|)) with the template-fitting photometric redshift Zfemp of Sawicki et. 
al. (1997)7for 848 galaxies in the HDF with / < 27 and measured UBVI. The 
symbols are the same as in Fig.l. The solid diagonal line indicates Za = ztemp', 
the dotted lines mark the region \za — ztemp\ < 0.5. About 90% of the galaxies 
fall within the dotted lines. 

There are 32 new spectroscopic redshifts of 2 < 4 HDF galaxies (Dickinson 
1998, Phillips et al. 1998, Steidel et al. 1999) listed in the recent paper by 
Fernandez-soto et al. (1999) which were not available to us when we derived 
our photometric estimator (Wang et al. 1998). Fig. 3(a) shows the difference 
between our estimated redshift Za and the spectroscopic redshift z, scaled by 
(1 -|- z), as a function of z. As can be seen, our empirical redshift estimator 
yields excellent predictions for z < 2, with a^ = 0.13 for all 18 z < 2 galaxies, 
and with cr^ = 0.08 if we omit the most discrepant object. For z > 2 galaxies, 
we find Gz = 0.34 for all 14 galaxies, and o"^ = 0.31 if we omit the galaxy with 
z = 2.002. The larger dispersion between Za and z for 2; > 2 is largely due to the 
small number of spectroscopic redshifts which were used in deriving the redshift 
estimator. 





Figure 3. (a) The difference between our estimated redsfiift Za (given 
by Eqs.(|q)-(|lO|)) and the spectroscopic redshift z, scaled by (1 + z), as a 
function of z, for 32 new z's. (b) Same as (a), but with the addition of 
Eq.(|ll|) to select flat spectrum galaxies with estimated redshifts given 
by Eq.(12). Star symbols represent the modified redshifts. 



3. Flat spectrum galaxies 

Galaxies in the redshift range \A <^ z <^ 2 have flat spectra, which makes the 
measurement of their spectroscopic redshifts difficult. The absence of calibrators 
in the range 1.4 < z < 2.2 (see Fig.l) leads to an artificial gap in Fig. 2(a) for 
1-5 ^ 2:a ^ 2 in our estimated redshift Za- 

To illustrate how the gap in Fig. 2(a) can be filled, we select galaxies which 
satisfy 



23.4 < / < 26.84, U > 26, 

\V -I\< 0.6, \B-V\< 0.8, 0.6 < [/ - B < 2. 



(11) 



We assume that these galaxies satisfy the same color redshift relation as the color 
range cr= 3 galaxies (which are closest to the intermediate redshift galaxies in 
color and redshift), but with a different constant offset (see Eq.(B)), i.e.. 



1.5 + 0.480(C/ -B)- 0.513(5 -V)- 0.250(F - /). 



(12) 



The constant offset value has been obtained by calibrating with a galaxy which 
satisfies Eq.(|^) and has spectroscopic redshift z = 2.002. Fig.2(b) presents 
our estimated redshift Za versus the Sawicki et al. template based photometric 
redshift ztemp] it is the same as Fig.2(a), but with the addition of Eq.(^) to select 
flat spectrum galaxies (represented by star symbols) with estimated redshifts 
given by Eq. ([T2| ) . This addition improves on the previous comparison by moving 
a number of galaxies with ztemp ^ 2 and Za ^ 2.5 or Za ^ 0.5 in Fig. 2(a) into 
the redshift gap of Fig. 2(b). 

We can modify our empirical photometric redshift estimator by adding 
Eq.([ll|) to select flat spectrum galaxies with estimated redshifts given by Eq.(12), 



before applying Eqs.(|6|)-(p^). Fig. 3(b) shows the difference between our mod- 
ified estimated redshift Za and the spectroscopic redshift z, scaled by (1 + z), 
as a function of z, for 32 newly measured galaxy redshifts. The modified red- 
shift estimator makes the estimated redshifts for z > 2 galaxies less biased, with 
Gz = 0.30 for 13 z > 2 galaxies (excluding the calibrator galaxy with z = 2.002). 

4. Discussion 

The advantage of our method is that it is model-independent and simple. The 
disadvantage of the method is that it requires a galaxy redshift sample for train- 
ing; this is clearly seen in Fig. 2(a), where an artificial gap for 1.5 ^ Za ^2 exists 
as a result of the absence of calibrators with measured redshifts in the range 
1.4 < z < 2.2 (see Fig.l). In this paper, we have illustrated how this gap in 
estimated redshifts can be filled by modifying the empirical photometric redshift 
estimator to include flat spectrum galaxies (see §3). It is important to fill the 
gap in measured spectroscopic redshifts for 1.4 ^ z ;^ 2, as it would enable 
better calibration of photometric redshift estimators, and better understanding 
of the nature of intermediate redshift galaxies. 

Since our method is purely empirical, it is sensitive to the data properties 
of the training set. In order to find and remove possible systematic effects 
in the photometric redshift estimator due to selection effects in the redshift 
measurements, it is important to have spectroscopic redshift surveys which are 
complete to a given magnitude limit. 

The photometric redshift estimator described here is simple and accurate. 
It can be further improved by additional data. 
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